Neutron Scattering Study on Spin Dynamics in Superconducting (Tl,Rb)2Fe4Se5 
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We observed in superconducting (Tl,Rb)2Fe4Se5 spin-wave branches that span an energy range 
from 6.5 to 209 meV. Spin dynamics are successfully described by a Heisenberg localized spin model 
whose dominant in-plane interactions include only the nearest ( J\ and J[) and next nearest neighbor 
(J2 and J2) exchange terms within and between the tetramer spin blocks, respectively. These 
experimentally determined exchange constants would crucially constrain the theoretical viewpoints 
on magnetism and superconductivity in the Fe-based materials. 

PACS numbers: 74.25.Ha, 74.70.-b, 78.70.Nx 



One astonishing aspect of the recently discovered T c ~ 
30 K iron selenide superconductors [l|-[5j is the coexis- 
tence of a large magnetic moment (3.3/ie/Fe) and high 
transition-temperature (T/v ~ 470-560 K) antiferromag- 
netic order 0, @|. Different from all other families of 
the Fe-based superconductors, the new iron selenide su- 
perconductors consist of Fe plates with highly ordered 
vacancy superstructure @-@|. Samples with less 
developed y/b x \/b order are known to not be supercon- 
ducting either when their compositions are close to Q or 
deviate significantly from [10] the ideal A2Fe4Se5 formu- 
las. In the latter case, an additional phase of orthorhom- 
bic vacancy order with a y/2 x y/2 or 2y/2 x a/2 unit 
cell also exists at intermediate temperatures. The large- 
moment block antiferromagnetic order developed on the 
vacancy ordered Fe lattice [Fig. l(a-b)] exists in super- 
conducting as well as insulating samples [Io| . 

The perfect >/E x >/E vacancy order demands one va- 
cancy per five Fe ions and the charge neutrality enforces 
a proportional number of intercalating A ions. While 
superconductivity can tolerate, or even requires, a small 
composition deviation from the ideal A2Fe4Se5 to dope 
charges, the crystal structure has to deform, as expected 
in any non-stoichiometric samples, such that the "va- 
cant" Fel site is found to be occupied by a few percent 
of Fe on average (g, 8] and fine-scale structural variation 
HI [T2| and phase-separation [13] are observed in super- 
conducting samples. Nonetheless, the block antiferro- 
magnetic order not only coexists with superconductivity 
in the same sample 0, 0, [HI EI ] , but also the interac- 
tion between these two long-range ordered states reveals 
itself in an anomalous magnetic order parameter near T c 
[si 0? [H| • Irrespective of the current theoretical debate 
on the role played by spin excitations in forming the su- 
perconducting Cooper pairs it is important to 
investigate spin dynamics experimentally in the new iron 



selenide superconductors to bound the discussion. 

Here we report inelastic neutron scattering measure- 
ments covering the whole spin excitation spectrum in the 
(Tl,Rb)2Fe4Se5 superconductor. Four doubly degener- 




FIG. 1. (color online) (a) Schematic diagram of the block an- 
tiferromagnetic structure in the I4/m unit cell. Only Fe ions 
with their spin directions are shown. J c is the exchange inter- 
action between spins in adjacent Fe planes, (b) Each shaded 
square highlights a block of four ferromagnetically coupled 
Fe 2 + ions in the Fe plate. The open (orange) and the crossed 
(blue) balls represent spins with opposite directions perpen- 
dicular to the afr-plane. The black line marks the unit cell. 
The four unique in-plane exchange interactions considered in 
this work are labeled, (c) Theoretical spin wave dispersions 
calculated using experimentally determined parameters. 
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FIG. 2. (color online) Constant energy slices of the acoustic branch of the spin wave excitations projected on the (if, K, 0) 
plane. The energy transfer is specified on each figure. The relative intensity is indicated by the color scale. The sample was 
aligned on one of the two crystalline twins in the 14/ m unit cell. The data were collected at SEQUOIA with Ei — 50 meV for 
(a-b) and 100 meV for the other panels. 



ate spin wave branches, one acoustic and three optical, 
form in three groups and span an energy range up to 
~210 meV [Fig. 1(c)]. A Heisenberg model, involving the 
intra- and inter-block nearest and next nearest neighbor 
interactions in the Fe plane and the nearest neighbor in- 
teraction between the planes, is sufficient to describe our 
data. Consistent with ab initio theoretical works, there is 
no need to resort to the third nearest neighbor exchange 
interaction J3. 

Single crystals of (Tl,Rb) 2 Fe 4 Se 5 (T c « 32 K) were 
grown using the Bridgman method [5|. A small single- 
crystal was used in a neutron diffraction study to de- 
termine the crystalline and magnetic structural proper- 
ties Q. For the inelastic studies, 240 plate-like crystals 
were mutually aligned on Al plates using an X-ray Laue 
diffractometer. The final assembly has a net sample mass 
of 19.5 g. The sample was sealed with He exchange gas 
inside an Al can [19| and measured with the SEQUOIA 



20|, [2l| fine resolution and the ARCS [22| wide angu- 



lar range Fermi chopper spectrometers at the Spallation 
Neutron Source (SNS) at Oak Ridge National Laboratory 
(ORNL). Neutron beams of incident energy Ei = 50, 100, 
200 and 350 meV were provided by the coarser resolution 
Fermi chopper [H| spinning at 180, 240, 360 and 420 Hz 
respectively on SEQUOIA. For ARCS E { = 400 meV was 
provided by the 700 meV Fermi chopper spinning at 420 
Hz. The sample was kept in its ground state by a closed 
cycle refrigerator operating at T w 6K. We will label the 
wavevector transfer Q = (iJ, K, L) using the tetragonal 
14/ m unit cell [6] of a = 8.683 and c = 14.39A. 



The importance of the nearest and next nearest neigh- 
bor exchange interactions in the Fe plane was identified in 
the initial ab initio study examining the iron pnictide su- 
perconductors . The lattice tetramerization forms the 
spin quartet block [Fig. 1(b)] making the intra- and inter- 
block exchange interactions inequivalent [24[. Therefore, 
the effective Heisenberg Hamiltonian 



(1) 



is used which includes five exchange constants Ji, J2, 
J{, J' 2 and J c as depicted in Fig. l(a-b), and the single- 
ion anisotropy constant A that quantifies the observed 
Fe spin S = 3.2(l)/g alignment along the c-axis 0]. 
This spin model on the \/b x Vb vacancy ordered lat- 
tice has been theoretically investigated (2^, [26[, and 
as to be shown later, describes the spin dynamics of 
(Tl,Rb) 2 Fe 4 Se 5 . 

Figure 2 shows the evolution of the acoustic branch 
of spin waves with increasing energy in (Tl,Rb)2Fe4Se5. 
The orientation of the tetramer Fe block with respect to 
the 14/ m unit cell can be clockwise [Fig. 1(b)] or counter- 
clockwise. The corresponding twins lead to eight, instead 
of four, Bragg spots for the magnetic {1,0,1} peaks pro- 
jecting on the (iJ, 0) plane. As energy transfer is in- 
creased, these spots broaden [Fig. 2 (a-b)] and then de- 
velop into well-resolved circular rings above ~ 30 meV 
[Fig. 2(c-g)]. The spin wave dispersions along the high 
symmetry direction [100] and [110] are demonstrated in 
Fig. 3(a) and (b), respectively. There exists featureless 
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nC 3. (color online) Slices of the spin wave spectrum (a) 
along the [100] and (b) the [110] direction, measured with 
Ei = 100 meV. The solid line is the theoretical spin wave 
dispersion described in the text. The spectral weight of the 
acoustic branch along the (c) [100] and (d) [110] direction was 
obtained from the constant-^ cuts at different L values. The 
solid lines are calculated intensities at corresponding L. 



scattering below 40 meV whose intensity increases with 
Q and is attributed to multi-phonon scattering. 

Across an energy gap, and above the acoustic branch, 
are two optic branches. Another gap proceeds a third 
optic branch at higher energy. The projection of the 
optical modes at 110±3 meV and 205±15 meV on the 
(if, if, 0) plane are shown in Fig. 4(a) and (c), respec- 
tively. To obtain eigenvalues of the optic modes at high 
symmetry points, constant- Q cuts at peak and back- 
ground positions were performed and their difference was 
fit to a Lorentzian. Example curves and fits are shown 
in Figs. 4(b) and (d). The peak at E = 209(1) meV in 
Fig. 4(d) presents the highest energy magnetic excitation 
mode in (Tl,Rb)2Fe4Ses. 

At the low energy end, the single-ion anisotropy A 
in Eq. (1) breaks the Heisenberg spin rotation symme- 
try, thus, opening a gap in acoustic spin waves at mag- 
netic Bragg points. The inter-plane coupling J c , which 
stabilizes the antiferromagnetic order at finite tempera- 
ture, also introduces a modulation in spin wave disper- 
sion along the c-axis. Fig. 5(a) and (b) show the details 
of the low-energy spin wave excitations obtained with the 
finer ^-resolution spectrometer configuration of Ei = 50 
meV. The energy gap in magnetic excitations is obvious. 
Constant- Q cut through the magnetic zone center (back- 
ground) is shown in Fig. 5(c). The difference intensity 
was fit to a step function convoluted with the instrument 
resolution to obtain the intrinsic gap value 6.5(3) meV. 
In Fig. 5(b), the dispersive curve of bandwidth ~18 meV 
along the c-axis is clearly observed. 

The simultaneous fit of the data from all branches 
along the multiple symmetry directions measured in this 
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FIG. 4. (color online) Constant energy slices of the op- 
tic branches of the spin wave excitations projected on the 
(if, K, 0) plane at the energy transfer of (a) E = 107 — 
113 meV and (c) E = 190-220 meV. (b) The constant-Q cut 
at (1.5,0.5,0) with a background at (3,-1,0) subtracted, (d) 
The constant-Q cut at (2,1,0) with a background at (3,1,0) 
subtracted. 



experiment to the spin wave solution of Eq. (1) yields the 
parameters: 



SJi = -30(l)meV, 
SJ 2 = 10(2)meV, 
SJ C = 0.8(l)meV, 



SJ[ =31(13)meV, 

SJ' 2 = 29(6)meV, 

SA = 0.3(l)meV. (2) 



The resulting spin wave dispersion curves in various high 
symmetry directions are shown in Fig. 1(c). They are also 
reproduced as the solid lines in Fig. 3(a-b) and Fig. 5(a- 
b), and are in excellent agreement with the measure- 
ments. To further check the reliability of the fits, the 
inelastic neutron scattering intensity was also calculated 
using these fitting parameters and over-plotted with the 
observed intensity in Figs. 3(c) and (d). The theory 
agrees well with experimental results. 

The Ji, J{, J2, J c and A in Eqs. (2) have the cor- 
rect sign to stabilize the observed block antiferromag- 
netic order, while the weaker antiferromagnetic J 2 frus- 
trates the ferromagnetically aligned spin block [Fig. 1(b)]. 
The strong difference of the exchange constants be- 
tween the intra and inter-block nearest neighbor and 
next nearest neighbor Fe spin pairs highlights the elec- 
tronic consequence of the lattice tetramerization in the 
y/h x \/b structure uncovered in structural refinement 
studies jf| B and emphasized by electronic structure cal- 
culations [24], I!?], . In particular, the recent ab initio 
linear response theory concludes that the significant in- 
plane exchange interactions include only Ji, J2, J[ and 
J2, whose calculated values [27j, remarkably, agree with 
our experimental results in Eqs. (2) qualitatively. 
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FIG. 5. (color online) Slices of the acoustic spin wave branch 
along (a) the [110] and (b) the [001] direction, (c) Constant- 
Q cuts through the magnetic Bragg point (2,1,1) and the 
background point (1.3,0.3,1). The difference intensity curve 
was fitted to the spin wave excitations convoluted with the 
instrumental resolution to obtain the energy gap at 6.5(3) 
meV. 



The block antiferromagnetic order not only exists in 
superconducting samples, but also in insulating samples 
with a less ordered y/b x y/b vacancy structure. Experi- 
mentally this type of order is observed even for 29% filling 
of the vacant Fe site [To| and ab initio calculations at 25% 
filling support the stability of the block antiferromagnetic 
order [27j |. Recently, the spin dynamics in an insulating 
Rbo.89Fei.58Se2 sample with the block antiferromagnetic 
order were investigated with inelastic neutron scattering 
[29j| . The overall energy scale of spin dynamics is sim- 
ilar to the system under study in this paper. However 
the published analysis of that data provides a somewhat 
anomalous result as compared to our result and the afore- 
mentioned ab initio works. Specifically their data could 
be least-squared fit only by a model that includes an 
inter-block J3 term in Eq. (1), where the preponderance 
of other results finds no need for this term. This discrep- 
ancy can be reconciled in two ways: One is to introduce 
a subtle physical effect in the ab initio studies that makes 
the insulating sample subtly different from the supercon- 
ducting sample. The other is to consider a small overlap 
of the observed excitation with the Q-resolution tail of 
its twin. We minimized this latter effect in our data by 
using the fine Q resolution of SEQUOIA for our acoustic 
mode measurements. 

In addition to the observed block antiferromagnetic or- 
der 0, 0, [Tq| , many other magnetic order configurations 
are possible with the Heisenberg model Eq. (1) with in- 
plane exchange interaction extending only to the next 



nearest neighbors on the y/b x y/b vacancy ordered lat- 
tice. Theoretical work has investigated the stability of 
these spin configuration s I24i426| . According to the cal- 
culated phase diagram [26|, the exchange constants de- 
termined in the present study put (Tl,Rb)2Fe4Se5 near 
the boundary between the block antiferromagnetic phase 
and a non-collinear antiferromagnetic phase. Further- 
more, such a non-collinear antiferromagnetic order has 
recently been observed in a spin-flop transition from the 
block antiferromagnetic order at 100 K in TlFei.6Se2 (or 
Tl2.5Fe4Se5) [30]. Therefore we postulate that composi- 
tion tuning has pushed Tl2Fe4Se5 just across this phase 
boundary. 

The success of the Heisenberg localized spin model in 
describing spin dynamics in the (Tl,Rb)2Fe4Se5 super- 
conductor may be attributed to the fact that the ob- 
served saturated magnetic moment 3.2(l)/ie/Fe is very 
large [7]. Therefore the system is close to the local spin 
limit. In the iron chalcogenide Fei+^Tei-^Se^ super- 
conductors without the intercalating layer and Fe vacan- 
cies, the spin excitations belong to the itinerant class as 
those in the antiferromagnet Cr which are described by 
more complex theories than linearized spin-wave theory 
3l|, I22I. Even-though, the binary iron chalcogenides are 



not simple antiferromagnetic metal like Cr: Their par- 
ent compounds have a magnetic wave vector that cannot 
be accounted for by Fermi surface nesting, and have the 
largest saturated magnetic moment (2/ie/Fe) among the 
previously discovered families of Fe-based superconduc- 
tors [33|. 

Moreover, a diffusive spin excitation component in the 
binary iron chalcogenides was recently resolved as origi- 
nating from interstitial Fe-induced short-range spin pla- 
quettes that contain the same four-spin blocks as found 
in the block antiferromagnetic order [34]. Such fluctu- 
ating spin quartets have been shown to contribute pro- 
nounced features in the spin dynamics in the parent com- 
pound FeTei.i [35|. The close link among these antiferro- 
magnetic states discovered experimentally in iron chalco- 
genides with or without vacancies, was anticipated in a 
theory including both itinerant and localized electronic 
states 

In summary, in this inelastic neutron scattering work, 
we contribute fresh insights to the understanding of iron 
chalcogenide superconductors by determining the spin 
Hamiltonian for a new block antiferromagnetic order in 
the (Tl,Rb)2Fe4Se5 superconductor. Our results agree 
with the majority of theoretical studies that state that 
the dominant exchange interactions extend only to the 
next nearest neighbor Fe pairs in the plane. The block 
antiferromagnetic order is frustrated only by the intra- 
block next nearest exchange J2, the weakest among the 
four dominant in-plane exchange interactions. Combin- 
ing our experimental exchange parameters with theoreti- 
cal calculations shows that this system is near the bound- 
ary of the block antiferromagnetism regime. A unified 
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theoretical framework for all observed types of magnetic 
order in the Fe-based superconductors should progress 
based on the experimental findings presented here. 
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